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Abstract. The Rubinstein-Duke or repton model is one of the simplest lattice model of reptation for the 
diffusion of a polymer in a gel or a melt. Recently, a slightly modified model with hardcore interactions 
between the reptons has been introduced. The curvilinear diffusion coefficients of both models are exactly 
determined for all chain lengths. The case of periodic boundary conditions is also considered. 
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1 Introduction 

The reptation of a polymer in an entangled melt was stud- 
ied long ago by De Gennes [I] who predicted the polymer 
length dependences of the curvilinear and self diffusion co- 
efficients as well as the viscosity and the relaxation time. 
Later, Rubinstein introduced a lattice model for the poly- 
mer motion incorporating most of De Gennes' ideas of rep- 
tation |2 . The so called repton model was then generalized 
by Duke to take into account the case of charged polymers 
during gel electrophoresis allowing the determination of 
the drift velocity This repton model seems particu- 
larly well adapted for DNA gel electrophoresis when the 
pore size is comparable to the persistence length |H5]. 

The theoretical prediction of the viscosity dependence 
with the polymer length is in apparent conflict with the ex- 
perimental observations | HI5ir7ll5| . This discrepancy is also 
observed in numerical simulations of the repton model 2 . 
Furthermore, a lot of interest in the calculation of the self- 
diffusion in the repton model focussed on the long polymer 
limit. The next to leading order term of the self-diffusion 
was long debated due to discrepancy between analytical 
and numerical results |Ulll(Jlllllll2lll3) . 

At the same time, the diffusion of a polymer chain in 
small channels attracts an increasing interest since it ap- 
plies to a great range of experimental situations. Brochard 
and De Gennes ^3] considered the case of a flexible poly- 
mer in a channel large compared to the monomer size but 
small compared to the polymer length whereas Odijk ^5] 
studied the case of stiff polymers with a persistence length 
larger than the channel width. The recent experimental 
access to nanometer scale channels allows to study the 
crossover behaviour between both regimes [T?)HTT| . The 
transport of long flexible polymer chains through Carbon 
nanotubes bring the interest to channels of width compa- 



rable to the monomer size allowing the determination of 
the curvilinear diffusion coefficients [TQ. Furthermore, the 
case of polymer diffusion in porous media with nanometer 
scale holes has been recently studied with a slightly 
modified repton model |2U|. This model presents identical 
dynamical rules for the curvilinear motion of the chain 
than the repton model and motivated the present interest 
on the analytical determination of the curvilinear diffusion 
coefficients as function of the chain length. 

The paper is organized as follows. In section 2 the rep- 
ton model is presented with the slight modifications intro- 
duced by Guidoni et al. [201 ■ in section 3, the main result 
of the paper, the exact calculation of the curvilinear diffu- 
sion coefficients, is presented and compared with previous 
results. Section 4 focuses on the repton model with peri- 
odic boundary conditions. Finally, in section 5, we discuss 
the next to leading order corrections of the exact curvilin- 
ear diffusion coefficients and we give some conclusions. 



2 Description of the model 

In the Rubinstein or repton model [2], the polymer chain 
contains N beads or reptons. A configuration C is charac- 
terized by N — 1 variables corresponding to the existence 
(n = 0) or not (rj = 1) of a stored length between the two 
reptons i and i + 1 along the chain. The polymer diffu- 
sion in a melt or a gel is obtained by the motion of stored 
lengths along the chain. Duke [3] generalized the model in- 
corporating the diffusion along a spatial direction subject 
to an electric field in order to model the electrophoresis 
of a polymer in a gel. The variables Ti are then replaced 
by N — 1 variables Uj with Uj = corresponding to the 
existence of a stored length and Oi = ±1 to the direction 
of the chain along the field between the reptons i and i + 1 
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Pa Pc Pc Pc Pb 

Fig. 1. A) A two dimensional representation of a polymer 
chain with 7 reptons and the underlaying cell structure in the 
repton model. B) The curvilinear representation of the same 
chain with horizontal lines linking two reptons in the same cell 
and with inclined lines for two reptons in different cells. The 
arrows represent the possible moves of the reptons with their 
respective rates (p a ,Pb or p c ). C) The representation of a poly- 
mer chain in the model with hardcore reptons. The existence 
of a gap between two reptons represents a hole. The possible 
moves of the reptons and their respective rates are represented 
by arrows. The identical curvilinear dynamics of both models 
is evident comparing the rates in B) and C). 



in the absence of stored length. In the simple Rubinstein 
model, the number of configurations is 2 N ~ 1 whereas it is 
increased to 3^"* in the Duke model. We are interested 
in the curvilinear diffusion coefficient of a neutral poly- 
mer or in absence of an electric field. In those cases, the 
direction of the chain between two reptons is not relevant. 
Thus, we restrict ourself to the Rubinstein model in the 
following. This reduces considerably the number of config- 
urations to consider without limitations on the generality 
of the results. 

In the Rubinstein model, three different moves exist 
(see Fignj: a) the motion of an end repton outside its 
cell when this cell is occupied by the neighbor repton, b) 
the motion of an end repton from its own cell to the cell 
of its neighbor if both cells differ, and c) the motion of 
an internal repton. Cell occupancies need to be consistent 
with this move which means that the repton leaves a cell 
occupied by a neighbor and reaches the cell of the other 
neighbor. The moves a) and b) correspond to the exit 
and entrance of a stored length in the chain whereas c) 
is the internal motion of a stored length. The respective 
rates of the moves are p a , pb and p c . In general, those 
rates are proportional to the number of accessible cells 



for the moving repton. Thus, in d dimensions and with 
a cubic lattice, pb — p c and p a — 2dpb- However, the 
existence of larger or smaller monomers at the two ends 
of the chain may either reduce or increase the ratio Pb/p c - 
Similarly, in case of a polymer embedded in a gel, the 
gel structure in which the polymer diffuses may affect the 
ratio Pa/pb- The later ratio is responsible as we shall see to 
the equilibrium properties and especially the curvilinear 
length of the polymer chain. The former ratio only affects 
the dynamics of the polymer. 

Guidoni et al. introduced recently a slightly mod- 
ified model with hardcore reptons in one dimension. The 
cells may only be occupied by a single repton and two 
neighbor reptons are either in neighbor cells (t.; = 0) or in 
next neighbor cells (ji = 1). The description of the chain 
in terms of the variables Tj is thus identical to the repton 
model. Tj = 1 corresponds to a hole between reptons and 
Tj = to the existence of a stored length. From the repton 
model, only the length of the chain is slightly modified to 
take into account the repton size. The polymer diffuses 
thanks to the motion of stored lengths or holes. Further- 
more, as seen on Fig. ^ the dynamical properties of the 
curvilinear motion are identical for both models. As a con- 
sequence, the curvilinear diffusion coefficients are equal. 



3 Exact curvilinear diffusion coefficients 

The curvilinear positions of the different reptons along 
the chain are defined by Sj = si + aJ2k<i T k with Si the 
position of the first repton and dTj the length between 
the two reptons i and i + 1. In case of hardcore reptons, 
an extra term (i — l)a should be added to take into ac- 
count the repton lengths. It is also interesting to define 
s c — ^j^Si/N the curvilinear position of the center of 
mass as well as s m — (s% + sjy)/2 the middle position 
of the polymer chain. This middle position differs from 
the center of mass one for a particular configuration C: 
5s = s c — s m = a — 2 i)ri/2N. This difference leads 

to an internal force that drives the polymer as will be 
shown later. 

Due to the N — 1 variables t, and their two possi- 
ble values, there exist 2 Ar ~ 1 internal configurations C = 
{ T i\i=i N-i f° r the polymer chain. We define the prob- 
ability P(C, s, t) for the configuration C to have a curvilin- 
ear center of mass position s c = s at time t. Due to the cor- 
relations between the internal configurations of the chain 
along the time, the evolution of the global variable s alone 
is not sufficient to determine the curvilinear diffusion co- 
efficient. The initial condition P(C,s,t = 0) = P eq {C)8{s) 
is assumed with 5(s) the usual delta function and P eq (C) 
the equilibrium probability of the configuration C. From 
this initial condition, the averaged curvilinear position 
over the ensemble of configurations remains zero at all 
times. Furthermore, the averaged squared curvilinear po- 
sition increases linearly with time proportionally to twice 
the curvilinear diffusion coefficient D C (N) of the polymer 
chain with N reptons. Note that we may expect the av- 
erage position for a single configuration C to shift a finite 
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i 


w(d — > C) 





VaT\ +Pb(i - n) 


1 to TV — 2 


Pciji + Ti+1 — 2t;T; + i) 


TV — 1 


p a T N -l +Pb(l ~ Tjv_l) 


i 


TO(C — > Ci) 





Pa(l — Tl) + PbTl 


1 to TV — 2 


Pa(ji + T i + 1 - 2TiT i + 1 ) 


TV- 1 


Pail - Tjv_i) +Pi,Tjv_i 



Table 1. Rates w(d — > C) and w(C — > d) in the repton model. 
The variables Ti corresponds to the configuration C. 



value from zero due to the initial difference Ss between the 
middle and the center of mass of the chain. The knowl- 
edge of this shift in position is necessary to determine the 
curvilinear diffusion coefficients. In fact, it reflects the cor- 
relations inside the chain between the configurations j20] ■ 
In order to define the dynamical rules, we introduce the 
configurations Ci for i — to TV — 1 which only differ from 
the configuration C by the fact that both variables Ti and 
Tj+i have changed to 1 — n and 1 — Tj+i. The changes from 
C to Ci are the only allowed moves in the repton model 
and it corresponds to a stored length passing through the 
repton i+1 (sec Fig^). In order for the move to exist the 
variables Tj and Tj+i must differ. With those constraints in 
mind, the evolution equation for the probability P(C, s, t) 
is: 

r) N ~ l 

-P(C,s,t) = -^2w(C^C i )P(C,s ) t) (1) 

i=0 
JV-1 

+ w ( c * -» C)P(Ci, s + As{C), t) 

i=0 

where the rates w are given in Tabled for open boundary 
conditions and As{C) = (t^i —Ti)a/N corresponds to the 
curvilinear motion of the center of mass. The variables 
To = 1 — Ti and "Hv = 1 — tjv-i have been introduced 
for consistency. The rate p a corresponds to the terminal 
repton exploring a new cell, the rate pb to the entrance of 
the terminal repton into the cell of its neighbor and the 
rate p c to the motion of stored length inside the chain (see 
Figd- Those rates are identical to those of the model of 
Guidoni et al. |2*fi| . 

The determination of the curvilinear diffusion coeffi- 
cients depends on three steps: the determination of i) the 
equilibrium probability P eq (C), ii) the long time limit of 
the averaged curvilinear position of a configuration C and 
iii) the long time limit of the derivative of the averaged 
squared curvilinear position (in the last case the average 
stands for the position as well as the configuration aver- 
ages). The different moments of the curvilinear positions 
of the center of mass are defined by: 

J s k P{C, s, t)ds = (s k )(C, t)P(C, t) (2) 

where P(C, t) is the probability to find the configuration C 
at time t. The time evolution for this probability is given 



by Eq. where P(C,s,t) is replaced by P(C,t). From our 
choice of the initial condition, P(C, t) = P eq {C) at all times 
and satisfies Eq.^with dP/dt = 0. The equilibrium prob- 
ability P eq (C) is the product of an identical probability 
P(t) for all variables r^: 

JV-l 

P eq (C) = [] P(n)- (3) 

The probability P(t) may be deduced from the evolution 
rates of the terminal reptons p a and p\> independently of 
p c . The particular configuration with all Tj = leads to 
p a P(0) = p b P(l). Thus, we deduce P(l) = p a /{p a +Pb) 
and P(0) = Pb/(p a + Pb)- The curvilinear length of the 
chain L = Sjy — s% = (TV — l)aP(l) where the overlinc 
stands for an average over the configurations with the 
equilibrium probability P eq (C). This result compares with 
the equilibrium length TVa + L in the model of Guidoni et 
al. [201 where the TVa difference comes from the length of 
the TV reptons. 

From Eq. ^ we deduce a differential equation for the 
average curvilinear position of the center of mass: 

SL(8)(C,t) = -Y,u(C->C i )(a)(C ) t) (4) 

i=0 

W-l p in \ 

+ £ wiC, - C) [{8)(d,t) - As(C)} 

where we used the property that / s k P(C, s + As, t)ds = 
{(s- As) k )(C,t)P(C,t) with P{C,t) = P eq {C). In the long 
time limit, the averaged curvilinear positions (s)(C, oo) 
saturate and the left hand part of Eq. ^vanishes. The same 
line of arguments on the averaged curvilinear squared po- 
sitions leads to the determination of the curvilin ear dif- 
fusion coefficient D C (N). The long time limit of d(s 2 )/dt 
equals: 

N-l 

2D C (N) = 2 Y w{C -*Ci)As{C){s)(C,oo) 

i=l 

2 N-l 
i=l 

After some transformations, the long time limit of the 
squared positions (s 2 )(C, oo) cancelled out. Consequently, 
D c depends on the rates w (Table^l and the average curvi- 
linear positions (s)(C, oo). For the later, we assume the 
following general expression: 

N-l 

(s)(C,oc)=aJ2f(i)T i (6) 
»=i 

where is a function to be determined. Consider the 
configuration with all variables r j; = 0, Eq. 0] leads to 
/(l) = —/(TV — 1). For the configurations Ck with t, = 
except T k = 1, if 2 < i < N-2, we obtain /(i+1) = 2f(i)- 
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f(i - 1) and, for i = 1, f{2) = (A + 1)/(1) + (A- l)/N 
with A = (p a +pb)/p c - Note the particular case A = 1 
or p a + Pb = Pc for which f(i) — is a trivial solution. 
This particular case corresponds to the lack of correlations 
as discussed by Guidoni et al. [201 ■ In the case A ^ 1, 
f(i) = A(A)i + B(A) with A{A) = -2B(A)/N and 
B(A) = (1 - Z\)/[(A r - 2) A + 2]. From this expression, 
it is possible to deduce (s)(C. oo) — 2B(A)Ss. It is easy 
to check that the expression for (s)(C, oo) given by Eq. [S] 
satisfies Eq. 0]for all configurations C. 

Two limiting cases are of interest. The fast internal mo- 
tion of reptons compared to the end-repton motion corre- 
sponds to A — > oo. In this case, (s)(C, oo) ~ —2 5s/ (N—2). 
In the case of slow internal motion of reptons, A — > and 
(s)(C, oo ) ~ (5s. It is interesting to notice the change of 
sign in the prefactor between those two limiting cases im- 
plying a shift (s)(C, oo) in different directions for the same 
difference 6s between the center of mass and middle posi- 
tion of the chain. 

The exact curvilinear diffusion coefficient D C (N) is de- 
duced from Eq. and: 

N-l 

Y, w(C -> d) = 2(N-2 + 2A)p c P(0)P(l) (7) 

i=0 

N-l 

£ w(C -» Ci)(ri+ X - n) = Pc (A - l)(n - rjv-i) (8) 



i=0 



(ri - rjv-i){s)(C, oo) = 2aP(0)P(l)/(l) (9) 



leading to: 
£> c (iV) = 



PaPb a 2 £> P(0)P(1) 

Pa+Pb (N-2)A + 2 N-2 + 2A- 1 



(10) 



with I?o = PcQ< 2 the individual diffusion coefficient of a 
free internal repton. This result is exact for all polymer 
lengths N > 2 and for all sets of rates p a , pb and p c . Note 
the linear dependence of D^T 1 with the number of reptons: 



D-\N) 



N -2 + 2A- 1 



(11) 



D P(0)P(1) 

FigEl illustrates this linear dependence for three dif- 
ferent sets of parameters. The equality between the rates 
p a and pb implies identical equilibrium properties (P(l) = 
P(0) = 1/2). The individual diffusion coefficient Dq = 1 
is imposed by setting a = 1 and p c = 1 leading to an iden- 
tical slope for all sets of rates considered. The sets only 
differ by the value of A, with A = 1/2, 1 and 2 from top 
to bottom, leading to a shift of the different lines on Fig|3 
The numerical results are obtained from Monte-Carlo sim- 
ulations. The average of the curvilinear squared position 
as function of time for 10 5 chains with initial conditions 
distributed with the equilibrium probability P eq {C) are 
determined for different chain lengths. Linear fits allow 
us to determine the curvilinear diffusion coefficients. The 
numerical errors are smaller than the symbols on FigEl 

It is possible to understand Eq.^Jas follows. There ex- 
ists N — 2 equivalent internal reptons and 2 extremal rep- 
tons. The inverse diffusion coefficient reflects the contribu- 
tion of the two kinds of reptons. The free internal reptons 




Fig. 2. Inverse of the curvilinear diffusion coefficient D c 
as function of the number of reptons iV for different sets of 
rates. The lines correspond to the exact analytical calculations 
and the symbols to numerical simulations with {p a ,Pb,Pc) = 
(1/4, 1/4, 1) circles, (1/2, 1/2, 1) squares and (1, 1, 1) diamonds. 



have a diffusion coefficient Dq. The probability for them 
to move inside the chain is P(0)P(1), the probability that 
there exists a unique stored length on the sides of the rep- 
ton. This leads to the contribution (N - 2)/D P(Q)P(l). 
The contribution of the two extremal reptons may be split 
in two parts: an extremal repton exploring a new cell (with 
a rate p a and a probability P(0) to have a stored length on 
its side) and an extremal repton moving into the cell of its 
neighbor (with a rate pb and a probability P(l) to have no 
stored length on its side). Both contributions are identi- 
cal (p a P(0) = pbP(l)) and add up to give the contribution 
2Z\- 1 / z) oP(0)P(l) to the inverse diffusion coefficient D^ 1 
in Eq. El 

In a d-dimensional square lattice, the choice p a = p c = 
1 and pb = 2d is customary and reflects the 2c? possible 
directions for an extremal repton to explore a new cell 
compared to the single cell possibility for the other moves. 
The curvilinear diffusion coefficient is then 



D C (N) = 



2da 2 



(2d+l)[N(2d+l) -4d}' 



(12) 



Exact results already obtained by Guidoni et al. [201 are 
recovered: the uncorrelated case A = 1 leads to D C (N) = 
D o P(0)P(l)/N and the particular case N = 3 to D c ('i) = 
DqP(0)P(1)A/ (A + 2). However, as can be seen from 
Eq. EH the second order in 1 /N differs from the one pro- 
posed in [2D1 d ue to the improper account of the correla- 
tions in their calculations. For example, the limit p c —> 
or A — y oo leads to the consistent limit of the curvilinear 
diffusion coefficient D C (N) ~ a 2 P(0)P{l)p c /(N - 2) -> 
in contrary to |2()j . In contrary, the fast internal motion 
of reptons (A — > or p c — > oo) leads to a finite D c — 
P(0)P(l)a 2 /2 independent of N. 
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Fig. 3. Inverse of the curvilinear diffusion coefficient D% b as 
function of the number of reptons N for different probabilities 
P(l). The lines correspond to the exact analytical calculations 
and the symbols to numerical simulations with P(l) = 1/2 
circles, 1/3 squares, 1/4 diamonds, 2/3 up triangles and 3/4 
left triangles. 



Fig. 4. Curvilinear diffusion coefficient D c as function of 
the number of reptons N in log-log scales for different sets of 
rates. The lines correspond to the exact analytical calculations 
and the symbols to numerical simulations with (p a ,Pb,Pc) = 
(1/4, 1/4, 1) circles, (1/2, 1/2, 1) squares and (1, 1, 1) diamonds. 



4 Periodic boundary conditions 

In the following, we consider the case of periodic boundary 
conditions studied by van Leeuwen and Kooiman [211122) . 
Those conditions correspond to introduce two new vari- 
ables "Hv and To with tn — To- The rates w(Ci — > C) and 
w(C — » Ci) for i = to N — 1 have the same expression 
given in Table \T\ior i — 1 to N — 2. The number of in- 
dependent variables t% in the periodic boundary case is 
increased to N compared to the N — 1 variables in the 
open boundary case leading to an increased number of 
configurations by a factor 2. 

In the case of periodic boundary conditions, the num- 
ber Ni of non-zero variables Tj or Nq = N — Ni of stored 
lengths is conserved by the dynamics. In this respect, 
the dynamics is non-ergodic. All configurations C with 
the same number of stored lengths have the same prob- 
ability P(Nq) and the equilibrium curvilinear length of 
the polymer with N stored lengths is L p b — sn — si = 
(N - l)an = (N - l)aNi/N since only N - 1 variables n 
are present in sn — s\. Furthermore, 

N—i 

J2 -> d) = 2 P MN - N 1 )/(N - 1) (13) 

4=0 

N-l 

^ w{C -> Ci)(r i+ i - n) = p c (r N - r ) = (14) 



i=0 



lead to: 



Dfi^Nt) 



a % p c N x {N-Ni) _ D P(0)P(l) 



N 2 (N - 1) 



N — 1 



(15) 



where we replaced No/N by the probability P(Q) to have 
a stored length and Ni/N by P(l) = 1 - P(0) in anal- 



ogy with the open boundary case. Note that the peri- 
odic and open boundary conditions present to the same 
leading behavior for the curvilinear diffusion coefficient 
D c - D o P(0)P(l)/N in the long chain limit. This result 
is also obtained for the self-diffusion coefficient but with 
a different length dependence |21U22| . 

FigCH illustrates the linear dependence of the inverse 
curvilinear diffusion coefficient for different probabilities 
P(l) as function of the number of reptons N. The numer- 
ical results are obtained from Monte-Carlo simulations by 
linear fits of the average of the curvilinear squared position 
as function of time for 10 chains with initial conditions 
comprising Ni = NP(1) variables Ti — 1. The numerical 
errors are smaller than the symbols. 



5 Discussion and conclusion 

In this paper, we have determined exactly the curvilinear 
diffusion coefficients D C (N) of the repton model. The cases 
of open and periodic boundary conditions were consid- 
ered. The inverse curvilinear diffusion coefficients present 
a linear behavior with the number of reptons in both sit- 
uations. The next to leading order term in the curvilinear 
diffusion coefficient D c shows interesting properties : 



Arp c (AQ 
A)P(0)P(1) 



1 



2(^-1) 
NA 



(16) 



This term is positive for A > 1 and negative for A < 1. 
The later case would correspond to a chain with large 
end-reptons limiting their motion compared to the inter- 
nal reptons. On FigQJ D C (N) as function of A" on a log- log 
scale is presented for three different values of A, from bot- 
tom to top A — 1/2, 1 and 2. The change of concavity of 
the curves is representative of the sign change in Eq. 1161 
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1/N 

Fig. 5. Effective exponent v c as function of 1/N for three 
different values of A. From bottom to top, A = 1/2, 1 and 2. 
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A possible consequence of the next to leading order 
term concerns the determination of the effective exponent 
obtain by a linear fit of the diffusion coefficient in a log-log 
scale. Let us define a finite size effective exponent v c (N) 
for the curvilinear diffusion coefficient D C (N) as follows: 

_ dlnD c _ 2(A - 1) 

~ dhiiV _ 1 + (N-2)A + 2 { ] 

Depending on A, this exponent is either larger or smaller 
than the expected one v c = 1 in the long chain limit 
except for the particular case A = 1 where v c (N) = 1 
(see Fig|SJ. The non- monotonous behaviour observed for 
the self-diffusion D s rjSJ is not present for the curvilin- 
ear diffusion D c . This difference is due to the absence of 
an anomalous behavior for D c whereas D s presents cor- 
rections of order 1/N 1 / 2 due to fluctuations in the chain 
length |2H]. 
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